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Ab initio wavefunction based methods are applied to the study of electron correlation effects on 
the band structure of oxide systems. We choose MgO as a prototype closed-shell ionic oxide. Our 
analysis is based on a local Hamiltonian approach and performed on finite fragments cut from the 
infinite solid. Localized Wannier functions and embedding potentials are obtained from prior peri- 
odic Hartree-Fock (HF) calculations. We investigate the role of various electron correlation effects 
in reducing the HF band gap and modifying the band widths. On-site and nearest-neighbor charge 
relaxation as well as long-range polarization effects are calculated. Whereas correlation effects are 
essential for computing accurate band gaps, we found that they produce smaller changes on the 
HF band widths, at least for this material. Surprisingly, a broadening effect is obtained for the O 
2p valence bands. The ab initio data are in good agreement with the energy gap and band width 
derived from thermorefiectance and x-ray photoemission experiments. The results show that the 
wavefunction based approach applied here allows for well controlled approximations and a transpar- 
ent identification of the microscopic processes which determine the electronic band structure. 
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I. INTRODUCTION 

The proper treatment of electron correlation effects in 
molecules and solids stands at the heart of modern elec- 
tronic structure theorvi 2 * 3 .^. In the study of molecular 
systems at least, wavefunction based quantum chemistry 
is known to provide a rigorous theoretical framework for 
addressing the electron correlation problem. The stan- 
dard quantum chemical methods^ make possible the 
construction of approximate wavefunctions at levels of 
increasing sophistication and accuracy and offer thus a 
systematic route to converged results. Advanced wave- 
function based calculations can be routinely performed 
nowadays for small and medium size molecules. However, 
algorithms able to treat electron correlation effects in pe- 
riodic systems are still at their infancy. The simplest cor- 
relation method is based on the M0ller-Plesset equations 
and second-order M0ller-Plesset (MP2) perturbational 
schemes for solids were actually implemented by several 



groups 



6,7,! 



',10,11.12,13,14 



As an extension of MP2, ab ini- 



tio many-body Green's function techniques were devel- 
oped too 15 i 16 ' 17 i 18 i 19 . In addition, investigations based on 
coupled-cluster (CC) theory were initiatecli 3 .! 2 ^ 2 ^ 2 ^ 3 -. 

The main point when applying quantum chemical 
methods to correlation calculations is to make use of the 
local character of the correlation hole which is surround- 
ing an electron. The latter optimizes the Coulomb re- 
pulsion between electrons and its accurate description is 
the essence of the correlation treatment. Starting point 
is a Hartree-Fock (HF) calculation on top of which the 
correlation calculations are implemented. For insulators 
and semiconductors, the development of efficient compu- 
tational tools is facilitated by the use of optimally local- 
ized Wannier functions. Several orbital localization pro- 
cedures were proposed in the context of periodic HF cal- 



culations. They normally rest on the a posteriori trans- 
formation of the crysta l Bloch orbitals, see for example 
Refs. HHiUll^Hi!. Alternative approaches were 
elaborated by Shukla et al2&£*££& and in the Toulouse 
group 11 ' 30 , where the self-consistent-field (SCF) equa- 
tions are solved directly in the Wannier representation. 

Local correlation methods such as the incremental 
scheme of Stol l 31 ' 32 and the local Hamiltonian formal- 
ism of Horsch et al~ and Grafenstein et a/. 34 ' 35 ' 36 
were applied before to the rigorous determination of 
both ground-statei 3 -^ 6 .! 2 !^ and excited-statei^ 3 -^^^ 
properties of infinite systems. Materials under inves- 
tigation were crystalline LiH, LiF, and LiCli^^ 6 -! 2 - 7 -^, 
diamond 38,40 , silicon 3 ^, beryllium 3 ^, carbon, and boron- 
nitrogen chains^ 3 -, and ircms-polyacetylene 37 . In the 
present paper we extend these studies to the case of ionic 
oxide compounds and choose MgO as a prototype and, 
at the same time, relatively simple insulating oxide. We 
compute correlation induced corrections to the valence 
and conduction HF energy bands and find that our final 
estimate for the fundamental gap is in rather good agree- 
ment with the experimental data. Our analysis provides 
also a clear picture of the major correlation effects that 
are responsible for the reduction of the HF band gap. 
The present study on magnesium oxide should open the 
way to similar calculations for more complex materials 
such as the 3c? transition-metal oxides. 



II. COMPUTATIONAL APPROACH 

Periodic HF calculations were carried out for crys- 
talline MgO with the crystal package 41 . Wannier 
functions associated with the HF valence and low-lying 
conduction bands were determined with the Wannier- 
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Boys localization module 24 of the same program. For 
conduction-band states, this may be a somewhat tedious 
task. However, localized Wannier orbitals can also be de- 
rived in more difficult cases with band disentanglement 
techniques 2 ^. 

For computing the correlation induced corrections to 
the valence and conduction bands we adopt the same 
quasiparticle picture and local Hamiltonian formalism as 

employed in Refs. mmmmi^IilliiO. physical 
reasoning is as follows. When an electron (or hole) is 

added to the TV-particle system, its surroundings relax 
and polarize due to the additional charge. This response 
of the system lowers the energy it takes to add the ex- 
tra electron (or hole). The particle plus the modified 
surroundings move together through the system in the 
form of a Bloch wave and define a quasiparticle. At the 
same time, some of the correlation contributions which 
are present in the ground-state of the iV-particle system 
are no longer operative in the (-W+1) [or (JV— 1)] -particle 
system, because some excitations are now blocked. This 
effect is called loss of ground-state correlation and must 
be accounted for too 33 ' 38 . A basic feature of our approach 
is that the correlation treatment is performed on a finite 
cluster C cut from the extended periodic system. With 
the cluster-in-solid embedding technique elaborated by 
Birkenheuer et aZi 38 ' 42 i 43 , the cluster C is divided into an 
active region Ca that supports the local occupied and 
virtual orbitals entering the post-HF calculations and a 
spatial "buffer" domain Cb including a number of atomic 
sites (and basis functions) whose role is to provide a good 
representation for the tails of the Wannier-like orbitals of 
the active region Ca ■ Thus C = Ca + Cb ■ All Wannier 
orbitals centered in Cb (and also the environment) will 
be held frozen in the correlation calculations. 

The orbital set associated with the finite-size clus- 
ter is obtained from the data supplied by the Wannier- 
Boys localization module of CRYSTAij2i. The interface 
program 4 ^ written for this purpose yields in addition an 
embedding potential corresponding to the frozen environ- 
ment, i.e., the surrounding HF electron sea. The detailed 
procedure for constructing the embedded cluster and the 
associated orbital set is described elsewher e - 42 ' 43 . Here 
we merely give an outline of this. The core, valence, and 
low-lying conduction-band cluster orbitals are generated 
by projecting the crystal Wannier functions onto the set 
of atomic basis functions attached to the cluster region 
C, 

K(R))= J2 l/3)^>K(R)>, (f) 

P,a£C 

where a and (3 are Gaussian basis functions centered 
within the region C, Sg a represents the inverse overlap 
matrix for the basis set attached to C, \w n (H)) is a Wan- 
nier orbital with index n and centered in the unit cell 
corresponding to the lattice vector R, and \w' n (R)) is 
its projected counterpart expressed exclusively in terms 
of basis functions centered in C. The \w' n (R)) func- 
tions are neither normalized nor are they orthogonal to 



each other because of the projection procedure mentioned 
above. Therefore they are group- wise orthonormalized in 
the following order: active core, active occupied, active 
low-lying conduction-band orbitals, buffer core, buffer oc- 
cupied, buffer low-lying conduction-band orbitals. This 
way the contamination of the most important types of or- 
bitals is minimized. Group-wise orthonormalization ac- 
tually means to use Schmidt orthogonalization for the 
inter-group orthogonalization and Lowdin orthonormal- 
ization inside the groups. This set of orthonormal or- 
bitals will be denoted by \w' n (R)). For the construction 
of the variational space to be used in the subsequent cor- 
relation calculations we follow the prescription suggested 
by Pulay and Saebo^. So-called virtual projected atomic 
orbitals (PAO's) are generated from the Gaussian basis 
functions associated with the active region Ca by pro- 
jecting out the occupied and the low-lying conduction- 
band orbitals \w' n (R)) via a Schmidt orthogonalization 
scheme. Thereafter, the PAO's are Lowdin orthonormal- 
ized among themselves to facilitate their subsequent use 
in the correlation calculations. 

By making the buffer region Cb sufficiently large the 
original Wannier orbitals of Ca are well represented. For 
example, for the Wannier orbitals corresponding to the 
lowest four conduction bands of MgO and centered at 
a Mg site, the six ligands surrounding the given Mg ion 
should be included in Ca because large contributions arise 
not only from the Mg 3s and 3p basis functions but also 
from the nearest-neighbor (NN) oxygen 2s, 3s and 2p,3p 
components. Farther neighbors need not be included in 
the central region Ca but may be put into the buffer 
zone Cb because the weight of the longer-range tails is 
small. Whereas a high-quality description can be easily 
achieved for the Wannier orbitals centered in the active 
region, the representation of the Wannier functions in the 
buffer zone is less accurate. The impact of this deficiency 
on the correlation calculations, however, is completely 
compensated by an appropriate choice of the embedding 
potential. Actually, the Gaussian orbital representation 
V^f b of the embedding potential is constructed by 

V™ h = F^ s -F{P c U, a,(3eC (2) 

where F crys is the self-consistent Fock operator from the 
periodic HF calculation and F[P C ] is the Fock operator 
associated with the density operator 

occ 

P C = 2]TK,)«| (3) 

arising from all occupied orbitals \w' v ) which enter the 
subsequent correlation calculations explicitly. This way, 
the correlation calculations are effectively performed in 
an infinite frozen HF environment. 

The data concerning the occupied and virtual orbitals 
of the cluster is transferred via the CRYSTAL-molpro 
interface^ 4 , to the quantum chemistry program mol- 
pro2£. The same holds for the matrix representation 
F"J S of the self-consitent Fock operator of the periodic 
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host system. The embedding potential itself is con- 
structed according to Eqs. ^ and §5§ using the MA- 
TROP module of the molpro program package. 

Local electron removal and electron addition one- 
particle configurations can be defined in terms of the set 
of occupied and virtual orbitals localized within the spa- 
tial domain Ca- 

Pr^I*} and K+ 1 )=4J1>), (4) 

where cr p<j and c^ Rqa are annihilation and creation op- 
erators for the valence and conduction band cr-spin or- 
bitals |iDp(R)) and \w' q (R)}, respectively, and |<fr) is the 
single-determinant ground-state wavefunction of the N- 
electron system. For clusters which are large enough, 
the Hartree-Fock valence and conduction energy bands 
of the periodic crystal can be recovered by diagonalizing 
k-depcndent matrices of the following form: 

HZW ^ e lkR (C> - £o HF K^> • (5) 

R 

Here stands for the reference unit cell and E^ F denotes 
the ground-state HF energy of the neutral iV-electron 
system. The diagonal terms R = 0, n' — n in this expres- 
sion are directly related to the on-site excitation energies 
within the Koopmans approximation^, i.e., ionization 
potentials 

IP p H » = (^I^K; 1 ) - < F = -e™ > (6) 

and electron affinities 

EA«» = «+ 1 |^l<+ 1 ) = -eg* , (7) 

with the latter being negative in the case of MgO. The 
off-diagonal terms are the hopping matrix elements in a 
tight-binding representation, 

ffiQl) = (*2g\H-E**\*ffi,). (8) 

Starting from the Hartree-Fock energy bands we in- 
clude next the effects of electron correlations. Within 
the quasiparticle approximation we may introduce for the 
(N+ l)-particle system a wavefunction of the following 
form: 

\9%£) = exp(S)<h I< ,\*), ( 9 ) 

where is the ground-state of the N electron system. 
The operator S does not commute with cj^ and de- 
scribes the correlation hole. For a more detailed discus- 
sion we refer to Ref. [l|. Alternatively, we may write 

I^r^t 1 ) = fic RgJ $ >> where I*) is the SCF ground-state 
and £1 acts like a wave- or Moeller operator—. Therefore, 
when correlations are taken into account, the analogue 
of ([5|) becomes 

H m Ak)=J2^ R «^\H-E a \^), (10) 

R 



where Eq is the energy of the iV-particle correlated 
ground-state |\&) (see also Refs. The cor- 

relation hole around an added electron (or hole) consists 
of a short-range and a long-range part. The short-range 
part originates from intra-atomic and short-range inter- 
atomic relaxation and polarization effects. We construct 
the short-range part of the correlation hole by separate 
orbital optimizations. Thereby the Wannier orbital to 
which the extra electron (hole) is attached is kept frozen 
and the changes within a finite region around it are deter- 
mined by performing an additional SCF calculatio n 38 ' 39 . 
This is sufficient if electron correlations are weak or mod- 
erate. If they are strong we would have to take into 
account the fact that changes in the nearby surround- 
ings are decreased when correlation effects among the 
electrons in that neighborhood are accounted for. The 
long-range part of the correlation hole consists of long 
range polarization of the environment. The effect of long- 
range polarization on the diagonal Hamiltonian matrix 
elements is estimated in this paper by applying the ap- 
proximation of a dielectric continuum. With the known 
dielectric constant of MgO, eo = 9.7, we can directly de- 
termine the polarization energy of a charge ±e outside a 
sphere of radius R. In addition, differential correlation ef- 
fects related to the existence of a different number of elec- 
trons in the system's ground-state and in the (iV±l) ex- 
cited states are investigated by subsequent configuration- 
interaction (CI) calculations. 

III. CORRELATION CORRECTIONS TO THE 
BAND STRUCTURE OF MGO 

Magnesium oxide is a prototype closed-shell ionic ma- 
terial that crystallizes under normal conditions in the 
rocksalt structure. It is extensively used in materials sci- 
ence as substrate for the epitaxial growth of films of other 
compounds. In recent years, it has attracted renewed in- 
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FIG. 1: Hartree-Fock band structure of bulk MgO. The core 
O Is and Mg Is — 2p bands are not shown in the figure. 



4 




FIG. 2: Plot of the Mg 3s-like conduction-band Wannier or- 
bital after projection onto a [MgigOi4] cluster. There is sub- 
stantial weight at the NN oxygen sites. 



terest because of its use as tunnel barrier in magnetic tun- 
nel junctions 47 . A different area where MgO attracted at- 
tention is optoelectronics. It turns out that when alloyed 
with ZnO, depending on the precise chemical composi- 
tion, the band gap of the system can be tuned along an 
interval ranging from 3.3 to 7.8 eV— . The latter number 
represents the fundamental gap of MgO 49 . 

The valence bands of MgO have oxygen 2p character, 
whereas the lowest conduction bands are mainly related 
to the Mg 3s and 3p orbitals, with some admixture from 
the O 2s, 3s and 2p,3p functions. The Hartree-Fock en- 
ergy bands are shown in Fig.l. We employed for our 
calculations the lattice constant reported by Sasaki et 
at, a = 4.217 A—, and Gaussian-type basis sets from the 
standard crystal library. Basis sets of triple-zeta qual- 
ity, Mg 8-511—, were used for the cations and triple-zeta 
basis sets suplemented with polarization functions, O 8- 
411*—, were applied for the more polarizable O ions. At 
the Hartree-Fock level and with this choice of the basis 
functions, the fundamental gap of the system is 16.20 
eV, i.e., 8.4 eV larger than the experimental value. As 
illustrated in Fig.l, there is a separation of about 16 eV 
between the O 2s and O 2p bands. There is also a clear 
separation between the low- lying Mg 3s, 3p complex and 
the other conduction bands. Density- functional calcula- 
tions within the local density approximation and using 
the same Gaussian basis sets as for the HF calculations 
give a gap of 5.0 eV between the valence and conduction 
bands, nearly 3 eV lower than observed in experiment. 

Projected Wannier orbitals associated with the Mg 3s, 
3p conduction-band complex are plotted in Fig. 2 and 
Fig. 3. We note that the Wannier-Boys localization mod- 
ule of the CRYSTAL program yields a set of Mg sp 3 hybrids 
for the lowest four conduction-band states. In order to 
arrive at a set of s and x, y, z -like functions, we applied 
to those projected hybrids a Pipek-Mezey localization 



procedure^, as implemented in the molpro package 4 ^ 
In this localization scheme, the number of atomic or- 
bitals contributing to a given molecular-like composite 
is minimized. Whereas the Mg 3s and 3p Wannier func- 
tions have substantial weight at the NN ligand sites, the 
valence-band O 2p (and 2s) Wannier orbitals are much 
more compact and contributions from neighboring ions 
are not visible in plots like those shown in Figs. 2 and 
3. In the calculations reported here the norm of the pro- 
jected active orbitals is typically larger than 99.0% and 
never below 98.5% of the original crystal Wannier func- 
tions. 



A. Correlation induced corrections to the diagonal 
matrix elements 

Relaxation and polarization effects in the immedi- 
ate neighborhood of an oxygen hole were computed by 
separate restricted open-shell HF (ROHF) calculations 
on [03gMg3o] clusters. The multiconfiguration MCSCF 
module of the MOLPRO package was employed for this 
purpose. The active region Ca of the [O3gMg 30 ] clus- 
ter contains a central 2p 5 (or 2s 1 ) 0~ ion and four of the 
twelve nearest 2 ~ ligands. These four oxygen neighbors 
are all chosen to be in the same plane and we denote them 
as O^.-.O 4 (see Fig. 4). In addition, we include in the 
cluster C the NN cations and all the O's in the next co- 
ordination shell of each of the O z xy sites, as sketched in 
Fig. 4. These additional Mg and O neighbors represent 
the so-called buffer region Cb and ensure an accurate de- 
scription of the tails of the orbitals centered in the active 
region Ca ■ 

When performing the ROHF calculations, the oxygen 
hole orbital is kept frozen^. We also freeze in our calcu- 
lations the core-like Is, 2s, and 2p shells of all Mg ions. 
Orbital relaxation effects arc listed in Table I for both 2s 




FIG. 3: Mg p-like conduction-band Wannier orbital after pro- 
jection onto a [MgigOi4] cluster. 
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and 2p oxygen holes. The on-site relaxation effect is quite 
large, more than 2 eV. There is also a substantial relax- 
ation/polarization effect associated with the first oxygen 
neighbors. For the four O l xy ligands included in the ac- 
tive region Ca of the cluster, this effect amounts to about 
0.41 eV in the presence of an O 2s hole and to 0.40 eV 
in the presence of a 2p hole. In Table I, we multiplied 
these numbers by three because there is a total of twelve 
oxygens in that coordination shell. We obtain thus a 
good estimate for the relaxation and polarization effects 
up to the nearest O neighbors. The fact that these relax- 
ation/polarization effects are additive was checked by ex- 
tra calculations with smaller, double-zeta basis sets^i on 
two different clusters: a [03gMg3o] cluster including only 
the four O xy sites in the active region and a [OssMgas] 
cluster where all twelve first oxygen neighbors were al- 
lowed to polarize. 

The four lowest-energy (iV+1) conduction-band states 
imply Mg 3s 1 and Sp 1 electron configurations. ROHF 
calculations were performed for such configurations on a 
[MgigOss] cluster with a [MgOg] kernel as active region, 
Ca- Beyond the [MgOg] kernel, this cluster incorporates 
again all Mg and O ions in the first two coordination 
shells of the active ligands. The on-site relaxation effects 
associated with the addition of an electron in a localized 
Mg 3s or 3p Wannier orbital are vanishingly small. The 
relaxation effects at the adjacent O sites induce energy 
shifts of 0.80-0.85 eV, see Table I. In these calculations 
the open-shell active orbitals (Mg 3s 1 or Sp 1 ) were again 
kept frozen^. We note that the energetic effect is nearly 
the same for the 3s 1 and ip 1 conduction-band states. At 
the scale of Fig.l at least, it induces an uniform down- 
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FIG. 4: Sketch of the [OagMgso] cluster employed for the 
calculation of short-range relaxation effects on the O valence- 
band states. The shortest line segments are Mg-O "bonds". 
The so-called active region Ca includes a central O ion and 
four nearest oxygen neighbors O l xy in the "horizontal" plane, 
see text. These active ions are shown as small black spheres. 



TABLE I: Correlation induced corrections to the diagonal 
Hamiltonian matrix elements for the valence-band O 2s, 2p 
and conduction-band Mg 3s, 3p states. All numbers are in 
eV. Negative corrections induce upwards shifts of the valence 
bands and shifts to lower energies for the conduction bands. 



AH nn (0) 


O 2s 


O 2p 


Mg 3s 


Mg 3p 


On-site orb. relaxation 


-2.64 


-2.04 






NN orb. relaxation 


-1.23 


-1.20 


-0.81 


-0.84 


Long-range polarization 


-1.80 


-1.80 


-2.25 


-2.25 


Total 


-5.67 


-5.04 


-3.06 


-3.09 



wards shift of the center of gravity of the 3s — 3p band 
complex. 

The data listed in Table I indicate that the on-site 
orbital relaxation and relaxation and polarization effects 
at the nearest oxygen sites in the presence of an extra 
electron or extra electron hole results in a reduction of the 
HF band gap by about 4.05 eV, that is, more than 45% of 
the difference between the HF and experimental values. 
Large corrections are also expected to arise from long- 
range polarization effects. The long-range polarization 
energy of a dielectric due to the presence of an extra 
charge ±e can be expressed as AE(oo) — AE(R) — C/R 1 , 
where 



C 



co 



1 



2e 



fill 



e is the static dielectric constant of the material, and R 
defines a sphere around the extra charge beyond which 
the dielectric response reaches its asymptotic value eo- 
The energy increment AE(R) denotes the relaxation and 
polarization energy up to the radius R around the added 
particle. The constant C can be obtained by choosing two 
different radii R\ and i?2 where the quantities AE(R\) 
and AE(R2) are computed, see for example Ref. l35l . 
However, we adopt here a simpler approach. We cal- 
culate the corrections due to long-range polarization by 
using the experimental value for the static dielectric con- 
stant, eo = 9.7. Since relaxation and polarization effects 
related to the nearest oxygen neighbors were already ac- 
counted for, both, for the valence-band hole states and 
the conduction-band electrons (see Table I) and since the 
core-like electrons of the Mg 2+ ions can be ignored in 
these calculations, we set R as the average of the radii of 
the first and second oxygen coordination shells around a 
localized 2p hole or 3s(3p) electron: R — (aV2/2 + a)/2 
for the O 2p valence-band states, where a = 4.217 A is 
the lattice constant^, and R = (a/2 + av / 3/2)/2 for the 
conduction-band states. The corrections to the diagonal 
matrix elements of the Hamiltonian arc then 



AH nn {0) 



ep-l e 2 
2e R 



(12) 



about -1.80 eV for the O 2s/2p bands and -2.25 eV for 
the Mg 3s/3p bands. These numbers are also included in 
Table I. 
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Before we discuss band narrowing (or broadening) due 
to correlation effects which also affect the band gap, we 
consider the loss of ground-state correlations. The latter 
leads again to a shift of the center of gravity of the bands. 
As pointed out in the previous section, some of the config- 
urations that are present in the TV-particle ground-state 
are blocked when an electron is added or removed. We 
investigated such correlations by CI calculations with sin- 
gle and double excitations (CISD) and discuss first dif- 
ferential effects for the N and (N—l) states. Since the 
oxygen valence-band Wannier orbitals are rather local- 
ized, we designed a cluster with a single O ion in the 
active region. Around this central O site we added one 
shell of Mg ions (6 Mg's) and two shells of anions (12+6 
O's) to build the buffer region Cb- In the CISD calcula- 
tions for the N and (N—l) configurations we correlate 
the 2s and 2p orbitals of the central O ion. Thereby 
the occupancy of the hole orbital is kept frozen in the 
calculations for the (N — l) states, which is referred to 
as the frozen hole approximatio n 38 ! 39 . Sets of separately 
optimized orbitals were used for the hole states of the 
(N — 1) -particle system, as discussed above. We found 
that for a 2p hole the correction to the on-site matrix 
element of the Hamiltonian is AH nn (0) = 0.85 eV, i.e., 
the 2p valence bands are downshifted by 0.85 eV. For the 
O 2s hole states, this correction amounts to 0.99 eV. One 
would expect similarly an upwards shift of the conduction 
bands. However, the situation is somewhat different here. 
When an extra electron is attached to the Mg 2+ ion, it 
polarizes the closed shells of the core. This is the domi- 
nant effect now because Mg 2+ has no valence electrons. 
We may employ for our analysis the high-quality results 
obtained for a free Mg ion by Doll et al££. The correc- 
tion to the ROHF Mg + ^ Mg 2+ ionization potential was 
found to be 0.27 eV in Ref. [53|. A similar differential cor- 
relation effect is occuring for the conduction-band states 
in bulk MgO. Therefore the conduction bands are shifted 
downwards instead of upwards and a partial cancellation 
between loss of ground-state correlation effects for the 
valence and conduction bands is taking place. The net 
result is a slight increase of the gap between the valence 
and conduction bands, in the range of 0.5 eV. Among 
the different contributions discussed here, this appears 
to produce the smallest corrections to the gap. More 
advanced calculations for studying such differential cor- 
relation effects are left for future work. 

To summarize the results listed in Table I, relaxation 
and polarization effects in bulk MgO are responsible for 
a reduction of the Hartree-Fock gap by 8.1 eV, which 
represents about 95% of the difference between the HF 
and experimental values, 16.2 and 7.8, respectively. Im- 
proved agreement is expected between our results and 
the experimental data when applying higher-quality ba- 
sis sets. We mention in this context that a reduction of 
3.8 eV is obtained for the HF gap when going from va- 
lence double-zeta to triple-zeta basis sets. Nevertheless, 
this large energy difference is mainly related to the very 
poor representation of the conduction-band states in the 



TABLE II: Nearest-neighbor (NN), Raw = (1, 1, 0)o/2, and 
next-nearest-neighbor (NNN), Rjvjvjv = (1,0, 0)a, hopping 
matrix elements for the conduction-band Mg 3s and 3p or- 
bitals, see text. Results of frozen-orbital CI (FO-CI) are listed 
in the second column; NOCI results in terms of separately 
optimized, relaxed orbitals (RO-NOCI) are given in the third 
column. All numbers are in eV. 



tnn 1 (R) 



FO-CI 



RO-NOCI 



ijvjv : 
3s — 3s 

3Px(y) — 3p x ( y ) 
3Px(y) — 3Py(x) 

3p z — 3p z 

tNNN ■ 

3s — 3s 
3p x - 3p x 

3Py{z) ~ 3p y ( z ) 



0.41 
0.66 
0.72 
0.13 

0.36 
0.77 
0.13 



0.42 
0.69 
0.77 
0.13 

0.37 
0.74 
0.12 



calculations with double-zeta basis functions and such ef- 
fects will be less substantial by further extension of the 
basis sets. We also keep in mind that differential correla- 
tion effects due to the existence of a different number of 
electrons in the system's ground-state and in the (N±l) 
excited states determine a small correction in the oppo- 
site direction, i.e., a slight increase of the fundamental 
gap. 



B. Off-diagonal matrix elements 

We discuss next the effect of correlations on the widths 
of the different bands. For that purpose the off-diagonal 
matrix elements of the effective Hamiltonian (Tl0|) have to 
be determined, i.e., the so-called hopping terms. At the 
Hartree-Fock level, these matrix elements are obtained 
by solving 2x2 secular equations where both wavefunc- 
tions are expressed in terms of localized HF orbitals, see 
Eqs. ((3]), (O, and ©. Relaxation and polarization ef- 
fects in the nearby surroundings of the added electron 
(or hole) are obtained by separate SCF optimizations 
for the (Ail) states. The separate optimization of the 
(Ail) wavefunctions leads to sets of non-orthogonal or- 
bitals. There will be thus both Hamiltonian and over- 
lap matrix elements between the extra electron (extra 
hole) wavefunctions ^a^a an d The calculation 

of such matrix elements has been recently implemented 
in molpro^. It is based on the transformation of the 
CI vectors to bi-orthogonal orbitals and follows an idea 
suggested by Malmqvist 55 . A similar approach was devel- 
oped by Broer et al. in Groningen^. Instead of putting 
these matrix elements directly into the eigenvalue equa- 
tions that determine the Bloch energies e n (k) we extract 
from these data effective hopping parameters associated 
with various pairs of orbitals. Comparison between such 
effective hopping terms and the HF off-diagonal Hamil- 
tonian matrix elements offers an insightful picture of how 
correlation effects modify the inter-site interactions and 
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is: 



FIG. 5: Sketch of the [Mg2803g] cluster employed for the 
calculation of the NN conduction-band hopppings. The active 
region Ca includes two NN Mg sites, small grey spheres, plus 
the bridging and apical ligands, large black spheres. 



consequently the widths of the bands. For energetically 
degenerate states, the effective hopping is defined as 



tnn' — {H nn ' S nn iH nn )/(l S nn i) 



(13) 



where H nn i and S nn i are the Hamiltonian and over- 
lap matrix elements between the (N ± 1) states n and 
n' . Since the separately optimized wavefunctions are ex- 
pressed in terms of sets of non-orthogonal orbitals, this 
type of secular problem is usually referred to as non- 
orthogonal CI (NOCI). In the case of mutually orthogo- 



nal states S n „i is zero and t r , 



R = H„ 



(R) 



Nearest-neighbor and next-nearest-neighbor (NNN) 
hopping matrix elements for the more diffuse conduction- 
band Mg 3s and 3p orbitals are listed in Table II. In the 
rocksalt structure, there are two O ions bridging two NN 
cations, whereas two NNN Mg's share a single oxygen. 
We designed a [Mg 2 s036] cluster for calculating the NN 
hoppings and a [Mg3g047] cluster for the NNN matrix 
elements. As active regions, we employed [Mg 2 Oe] and 
[Mg20g] kernels, respectively, see Fig. 5 and Fig. 6. As 
in the calculations for the on-site matrix elements, for 
each of these clusters we included in the buffer region Cb 
all metal and O ions in the first two coordination shells 
of the active oxygens. Since the second cluster consists 
of 83 atoms, we reduced the computational effort by re- 
moving the polarization functions at the oxygen sites for 
this particular cluster. Two different values are given in 
Table II for each hopping matrix element. Results ex- 
tracted from 2 x 2 CI calculations in terms of frozen HF 
Wannier orbitals are listed in the second column. In the 
third column, we allowed for full relaxation of the 2s and 
2p orbitals of the O ions included in the active region Ca, 
for each of the (JV + 1) configuration state functions en- 
tering the 2x2 CI. As illustrated in Figs. 5 and 6, those 
active anions are the ligand(s) bridging two Mg sites and 
the ligands which are nearest oxygen neighbors of one of 



the active Mg sites and also of the bridging ligand(s). 

The results show that the separate optimization of the 
(JV+1) wavefunctions induces only minor changes on the 
electron hoppings, in the range of few percent. An inter- 
esting feature is that variations occur in both directions, 
i.e., some of the effective hoppings are slightly enlarged 
by taking into account short-range relaxation and polar- 
ization effects and some are reduced. Since these changes 
are quite small, the width (and the overall structure) of 
the lower conduction bands will change very little. Re- 
garding the longer-range polarization effects, we expect 
that their influence on the hopping terms is negligible. 

In the case of non-equivalent orbitals, neighboring Mg 
3s and Mg 3p, it is more difficult to define some NOCI ef- 
fective hoppings and the corresponding data is not shown 
in Table II. Nevertheless, the effect of short-range relax- 
ation and polarization is also small for these interactions, 
with reductions of the CI splittings of few meV. 

The same type of analysis was applied to the (JV— 1) O 
2s 1 /2p 5 hole states. The spatial extent of the oxygen or- 
bitals and the inter-site matrix elements are significantly 
smaller. Since the largest relaxation effects concern or- 
bitals at the same site, see Table I, we included in a first 
step only two ligands in the active regions of our clusters. 
Nevertheless, few other atomic shells were added around 
these active sites in the buffer region C b ■ We considered 
a [O28Mgi ] cluster for computing the relaxation effects 
on the NN hoppings and a [C^Mgn] cluster for the NNN 
terms. The results are collected in Table III. The largest 
corrections arise for the NN matrix elements, with abso- 
lute values that are similar to the corrections obtained for 
the NN conduction-band hoppings. In relative numbers, 
these corrections are somewhat larger, by 10% to 20% 
for the NN hoppings. A very different situation occurs in 
strongly correlated oxides such as the layered cuprates, 
where the existence of an antiferromagnetic spin back- 
ground determines a reduction of the effective quasipar- 
ticle hoppings by a factor of four—. 



=!i 



HI 



FIG. 6: Sketch of the [Mg3e047] cluster employed for the cal- 
culation of the NNN conduction-band hopppings. The active 
region Ca includes two Mg sites, small grey spheres, and nine 
O neighbors, large black spheres. 
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TABLE III: NN, Rjvjv = (l,l,0)o/2, and NNN, Rjvivjv = 
(1,0, 0)a, hopping matrix elements for the valence-band 
states. Results of frozen-orbital CI (FO-CI) are listed in 
the second column; NOCI results in terms of separately op- 
timized, relaxed orbitals (RO-NOCI) are given in the third 
column. All numbers are in eV. Data for the oxygen 2s or- 
bitals are also included, although the O 2s band is much be- 
low the O 2p bands. The values in parentheses include orbital 
relaxation effects at four additional O sites, see text. 



tjin' (R) 



FO-CI 



RO-NOCI 



tjvjv : 
2s - 2s 

2Px(y) ~ 1Py(x) 

2p z - 2p z 

tNNN '■ 

2s - 2s 
2p x - 2p x 



0.10 
0.32 
0.42 
0.12 

0.01 
0.06 



0.12 
0.37 
0.49 
0.14 

0.01 
0.06 



(0.11) 
(0.36) 
(0.47) 
(0.13) 



The trends displayed in Table III are confirmed by 
multi-reference CISD [MRCI(SD)] calculations. We cor- 
related in these calculations the eight 2s and 2p orbitals 
at the two NN oxygen sites and the reference active space 
included those two orbitals involved in the hopping pro- 
cess. The hopping matrix element is half of the en- 
ergy separation between the lowest two eigenstates. The 
MRCI values for the 2p x — 2p x and 2p x — 2p y hoppings, for 
example, are 0.33 and 0.43 eV. These numbers are again 
larger than the HF results, although the magnitude of the 
effect is smaller as compared to the NOCI calculations. 
Having the data of the SCF calculations for the (N— 1) 
oxygen hole states at hand, we performed an analysis of 
the changes produced in the composition of the (relaxed) 
orbitals in the immediate vicinity of an 2s 1 or 2p 5 
anion. We found that an oxygen hole is causing polariza- 
tion and "bending" of the 2p orbitals at the NN ligand 
sites. The "bending" of the NN 2p orbitals towards the 
0~ site takes place through both p—s and Pi~Pj mixing. 
These effects result in stronger inter-site orbital overlap 
and explain the fact that the effective NOCI hoppings 
are larger than the corresponding HF values. A similar 
analysis for the conduction-band (N+l) states is compli- 
cated by the presence of several sets of "active" orbitals, 
Mg 3s, 3p and bridging O 2p, and we could draw no clear 
conclusions in that case. 

Extra calculations were performed for the NN hop- 
ping matrix elements with four additional oxygen sites 
included in the active region of the cluster, Ca- Those 
are the four ligands which are nearest neighbors to both 
O ions involved in the hopping process. They are situated 
in the median plane of the segment R_/vat = (1, 1, 0)a/2. 
A cluster composed of 62 sites, [C^oMgaa], was employed 
for these calculations. The results are given in Table III 
in parentheses. The corrections due to relaxation and po- 
larization at the four nearest O neighbors are very small, 
0.01 to 0.02 eV. The effect of these corrections is to re- 
duce somewhat the absolute values of the NN hoppings. 



Experimental studies for characterizing the valence 
electronic structure of MgO have been carried out 
using x-ray photoelectron spectroscopy (XPS)££ and 
angle-resolved ultraviolet photoelectron spectroscopy 
(ARUPS)^. The measured width of the O 2p bands 
is about 6.5 e V 58 i 59 . The HF valence-band width is 5.50 
eV and inclusion of local correlations leads to a slight 
broadening of the O 2p bands, which brings the ab initio 
result in good agreement with the experiment. For com- 
parison, density-functional calculations within the local 
density approximation and using the same Gaussian ba- 
sis sets as in the HF calculations predict a width of 4.68 
eV for the O 2p bands. 

The correlation induced corrections to the widths of 
the bands also influence the band gap. In the fee lat- 
tice the dispersion of p bands at the T point depends on 
two of the nearest-neighbor hoppings, e x (T) = const. + 



8/ 



(110) 



4i 



(Oil) 



60.61 



ble III, fi ^ 



(110) 



and 0.01 eV for t x 01 x 1] 



With the notations from Ta- 
Corrections of 0.04 eV for t x \ x °^ 

(or e io) ), see Table III, imply an 



upwards shift of the O 2p bands at the T point and a 
narrowing of the fundamental gap by about 0.30 eV. For 
the Mg 3s— 3p conduction-band complex such changes at 
the T point are smaller because the correlation induced 
corrections to the 3s — 3s and 3s — 3p inter-site matrix 
elements are lower. 



IV. SUMMARY AND CONCLUSIONS 

We have analyzed the different correlation contribu- 
tions to the energy gap of MgO and to the widths of 
the conduction and valence bands. This was done within 
the quasiparticle description. As regards correlation ef- 
fects we have distinguished between relaxation and po- 
larization around an electron or hole added to the (neu- 
tral) ground-state. The net result is a reduction of the 
Hartree-Fock gap from 16.2 eV to a value of 8.1 eV. This 
has to be compared with a measured energy gap of 7.8 
eV. Within the local density approximation (LDA) to 
density functional theory a gap of 5.0 eV is found. This 
is not surprising since LDA is known to produce too small 
gaps for insulators. 

The calculations were performed with triple-zeta basis 
sets. Since on the Hartree-Fock level the calculated gap 
differs for double- and triple-zeta basis sets by 3.8 eV, one 
may consider the good agreement with the experimental 
gap as somewhat fortuitous. Triple-zeta basis sets are 
known, however, to produce reliable results in quantum 
chemistry and therefore a further extension of the basis 
set should keep the corrections small. Presently we are 
not able to work with larger basis sets. 

It was shown that a large contribution to the correla- 
tion induced corrections to the fundamental gap comes 
from on-site and nearest-neighbor relaxation, i.e., from 
the immediate neighborhood of the added particle. But 
also the long-range part of the polarization generated by 
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the extra particle contributes substantially to the reduc- 
tion of the HF gap. This long-range part can be treated 
in a continuum approximation thereby using the known 
dielectric constant of MgO. The so-called loss of ground- 
state correlations makes a small contribution in MgO. 
The reason is that the conduction-band Wannier orbitals 
have predominant Mg 3s or 3p character. The added 
electron will essentially go thus to a Mg site where it 
polarizes the closed Is 2 , 2s 2 , and 2p 6 shells. This effect 
reduces the gap and counteracts the loss of ground-state 
correlations which occurs when an electron is removed 
(hole state) and therefore can no longer contribute to the 
correlations of the remaining ones. Finally, also changes 
in the widths of the bands influence the energy gap. 

One surprising effect which we found is an enhance- 
ment of the width of the valence bands when local corre- 



lations are taken into account. This is the opposite one 
expects since correlations lead usually to band narrowing 
instead of broadening. It results from a slight bending 
of nearest-neighbor 2p orbitals towards a 0~ site when 
we account for correlations. The bending of the nearest- 
neighbor ligand orbitals increases the wavefunction over- 
lap and hence facilitates the hopping. 

The present results prove the usefulness of band cal- 
culations based on quantum chemical techniques. They 
allow for well controlled approximations and a transpar- 
ent interpretation of the different microscopic processes 
which determine the size of the energy gap and the widths 
of the bands. 

We thank Dr. E. Pahl and Prof. J. Fink for fruitful 
discussions. 
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